tt a tt St Ht tt HH HH HH Ha 
tt a tt a a tt i 


T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 


HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HH HHH 
HHHHHHHHHHHHHHH 
HHHHHHHHHHHHHHH 
HHHHHHHHHHHHHHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 


tt tt tt st 4 Ht Ht Ht Ht Ht Ht HH HH 
AAA a 
| LE) 


LLLLLLLLLLLLLLL 
LLLLLLLLLLLLLLL 
LLLLLLLLLLLLLLL 


**F 1LE**]D**MTHDCOSH 


MM MM TTTTTTTTTT HH HH DDDDDDDD cccccccc 000000 SSSSSSSS HH HH 
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-TITLE MTHSDCOSH 3 Sovete Fleeting Hyperbolic Cosine routine 
- IDENT /1-008/ : File: MTHDCOSH.MAR Edit: RNH1008 
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COPYRIGHT (c) 1978, 1980, 1982, 1984 B 
DIGITAL EQUIPMENT CORPORATION, MAYNARD. MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


® 
g 
® 
® 
® 
THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY SE USED AND COPIED + 
ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
INCLUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER * 
COPIES THEREOF NAY NOT BE PROVIDED_-OR OTHERWISE MADE AVAILABLE TO ANY * 
OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY * 
TRANSFERRED. * 
® 

THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
AND SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT * 
CORPORATION. * 
+ 
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® 

® 
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DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 
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FACILITY: MATH LIBRARY 

“ABSTRACT: 

MTHSDCOSH is a function which returns the double Sloat ing mygornet is cosine 
s 


of its double precision floating point argument. The call standard 
call-by-reference. 


VERSION: 01 
HISTORY: 
AUTHOR: 
Peter Yuo, 29-Jun-77: Version 01 


MODIFIED BY: 
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MTHSDCOSH ; Double hee ting Hyperbolic Cosine rout 16-SEP-1 1:15: AX/VMS Macro V04-00 
1-008 HISTORY ; Detaile LS tetrent bale History 6=SEP-1 1382 er 221: 33 EMTHRIL. SR SRE SMTHDCOSH.MAR: 1 
3 6 9 -SBTTL HISTORY ; Detailed Current Edit History 
8 : g ; ALGORITHMIC DIFFERENCES FROM FP-11/C ROUTINE: none 
4 5 ; Edit History for Version 01 of MTHSDCOSH 
000 5 3 0-2 pTupcennon ,Shenges % to sprees ew . 
000 58, MTH$_... change 
44 59 ; Changed error A "eachinten. Put error result in RO:R1 before 
4 ¢ 3 calling MTHS$SSIGNAL in order to allow user modify error result. 
$009 6¢ 30-3 Seven term Voyt or series, in powers of orquaset:_t replaced 
0000 65 ; uy six term Chebyshev series, in powers o anger? 
0000 64 ; th overhang, to improve accuracy. 18-May-1978; Mary Payne 
0000 65 ; 0-4 Fix LOG(2) constant. TNH 16-June-78 
0000 66 ; 1-001 = Update version number and copyright notice. JBS_ 16-NOV-78 
0000 67 ; 1-306 - ¢ ange MTH_ FLOOVEMAT to MTHSK_FLOOVEMAT. JBS 07-DEC-78 
0000 68 ; 1-005 - Removed SSRADEF - not needed. “JBS 16-DEC-78 
0000 69 ; 1-004 = Add to the PSECT it tive. JBS 22-DEC-78 
0000 70 ; 1-005 = Declare externals. *,! May-1979 
0000 71 ; 1-006 - Use MTHSDEXP_R6. SBL 7-Sept- -1979 
4464 i 3: 1-007 - rer ed lower Limit for Chebyshev approximation from 2**-27 to 
0000 74; - Eliginated second call to exe for he yolues Pe oe -25 and 
0000 2 3 8.5*ln2 by gee COSH(x) = (2 + 1/2)/2, with Z = EXP(ix!). 
0000 76 ; - Eliminated second call to EXP for a. values Seeueen 28.5*in2 
0000 3 and 127*\n2. 
0000 78 ; - Changed all final floating point divisions by 2 to interger 
0000 79 ; subrtracts of 1 from the exponent field. 
0000 80 ; - Changed onery mask to excluded R7 - no Longer needed. 
0000 B) 3 - Extended maximum range from 87.69 to 128*\n2=88.72. 
0000 82 ; - phaneged sogte for computing EXP(ixi-ln2) to reduce error. 
0000 83 ; - RNH 10-FEB- 
0000 84 ; 1-003 > Changed W* to G* on call to MTHS$SIGNAL RNH 09-Sept-1981 
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; Declarative Part of Modul 6-SEP- 1:21:5 MTHRTL.SRCJMTHDCOSH.MAR; 1 


om 
r Leow 
“o 


I 
-SBTTL DECLARATIONS ; Declarative Part of Module 


INCLUDE FILES: 


t EXTERNAL SYMBOLS: 


. GBL 

-EXTRN MTHSDEXP_R6 
-EXTRN MTHSK_FLOOVEMAT 
-EXTRN MTHSSSIGNAL 


SOOQOCOoCoOooooooooooooooooooo }-— 
SOOSCOCOOCOCOCOOOCOOCOOSOOOCSOOSOoOSoOSO -8 
OOCooooCoooooooooCoooooooooo 


| 
| 
| 
| 
: EQUATED SYMBOLS: 


8 
9 
0 
4. 
38 
94 
95 
96 
97 
98 
99 
100 
101 
106 
00004080 10 $0.1.0 =z “F1.0 3; 1.0 
00000004 He value =4 3; value.rd.r 
106 ; 
107 ; MACROS: none 
108 ; 
109 ; 
4 449 : PSECT DECLARATIONS: 
00000000 \N¢ -PSECT _MTHSCODE PIC,SHR,LONG, EXE ,NOWRT 
0000 =—11 3 program section for math routines 
0000 114; 
0000 115 ; OWN STORAGE: none 
0000 116; 
0000 117; 
0000 118 ; CONSTANTS: 
0000 119; 
0000 120 
0000 121 D_127_L0G_2: 
2BDAC7E2 OF 334380 0000 1$¢ -GQUAD “X2BDAC7E20F 334380 : 127*Ln2 
0008 \$ D_128_L0G_2: 
CF79F7D1 72174381 Baye 124 -GQUAD  *XCF79F 70172174381 > 128*Ln2 
001 125 D_2_POWER 28.5: 
DE6433F9 O4F34EBS 001 126 J *XDE6433F 904F 34EB5 > 2**28.5 
001 127 D_LOG_2_HI: 
CF80F7D1 72174031 49 ' $ ry . “XCF BOF 70172174031 : (high 49 bits of Ln2) + 2%#-49 
: } 
FF81898C B6C3ASCA 09 180 “= "QUAD “XFF8189B8CB6C3A5CA ; ln2 = D_LOG_2_HI1 | 
0 : 1 § DCOSHTAB: | 
2853 3594 00 1 -WORD “00 $ 8h «0088183 
DBD2 EC 02C 134 WORD *00 4,°01557 3 DECIMAL: 0.2759648863787355D-06 
OCF 1 0° 030 135 -WORD *0034 $8 » 000639! 
CA4F AB2 034 #1 § wORD “01 2? 0145117 : DECIMAL: 0.2480155975461668D-04 
0860 6 0038 1 WORD *00 §.20003340 
FEO 8 1 Q003c 138 WORD *0134001,°01 8 ; DECIMAL: 0.1388888889781712D-02 | 
AAAA 3E2A v4 139 WORD *0037052,°012525 
6FC4 AAAA 44 140 WORD %*0125252,°0067704 3; DECIMAL: 0.4166666666665359D-01 
eit 4000 0048 141 WORD *%*0040 3 9 
005 0000 004C 142 wOoRD 0,*000006 ; DECIMAL: 0.5000000000000000D0 
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ering tecten Gets acaut EREEIONS CUDE:EY RTE MESRABEGRR ass PO 
3 WORD “yo ,0 
5 DCOSHLEN- = DCOSHTAB/8 


10N 


; DECIMAL: 0.1000000000000000D+01 
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p++ 
: FUNCTIONAL DESCRIPTION: 


1:33:38 pals Macro V04-00 


- Standard Double Precision Floating DCOSH 


MTHRTL.SRCJMTHDCOSH.MAR; 1 


: DCOSH = double precision floating point function 
DCOSH(X) is computed as: 
If ixi < 2e#-28, DCOSH(X) = 1. 
IF 2e*-28 =< {Xi < 9-23 COSH(X) = Chebyshev series 
If 0.25 =< ixi =< 28. «(n2, set Z = DEXP(:X!) and compute 
DCOSH(X) = (2 + 1/2)/2, 
If 28.5*ln2 < {Ki =< 127#ln2 COSH(X) = EXPCIK!)/2 
If seeing < iKi =< 128 ln2 then DCOSH(X) = DEXP(!X!-LOG(2)). 
If 128*ln2 < ix! then overflow. 
CALLING SEQUENCE: 
DCOSH.wd.v = MTHSDCOSH(x.rd.r) 
INPUT PARAMETERS: 
LONG = 4 3; define rt multiplier 
x = 1 * LONG ; Contents of x is the argument 
IMPLICIT INPUTS: none 


OUTPUT PARAMETERS: 
VALUE : 
IMPLICIT OUTPUTS: 
COMPLETION CODES: 
SIDE EFFECTS: 
Signal: MTHS FLOOVEMAT if 128*Ln2 < {Xi 


none 


none 


double precision floating hyperbolic cosine of the argument 


with reserved operand in RO/R1 


; (copied to the signal mechanism vector CHFSL_MCH RO/R1 Y LIBSSIGNAL). 
LIBRARY" 


; Associated message is: 


“FLOATING OVERFLOW IN MATA 


; operand -0.0 unless a user supplied (or any) error handler changes CHFS$L 


; NOTE: This procedure disables floating point underflow, enables integer 


low. 


-ENTRY MTHSDCOSH, “M<IV, R2, R3, 


MTHSFLAG_JACKET : 
G*MTHSSJACKET_HND, (FP) 


MOVAB 


R4, RS, R6 > 
standard call-by-reference entry 
disable DV (and FU), enable IV 


flag that this is a jacket procedure in 


set handler address to jacket 


. Result is reser 


handler 


1 3 | 
MTHSDCOSH ; Double Floating Hyperbolic Cosine rout 16-SEP-1984 01:15:5 AX/VMS Macro v04-00 P 
' WTHEDCOSH. \o Stotdard Double Poacteion es ecsepedobe Giiskie? PATANTy emeserviococu.man:1 2 6, 


; case of an error in routine 


1 
1 
1 ¢ : 
3 ; If an error, convert signal to user PC 
; ang resigna | 
$0 06 BC 70 Oos1 508 ovo (AP), RO aoe alee Tat = avelue(AP) | 
value 3 = iX!t = a@value 
0 4000 Br mA 0065 206 Bitu2  #xKBOG0" RO : ROURI = IX! 
FBO BF 0 =B1 6A Pw RO #°x3F 80 3 compare {Xi with 0.25 
14 18 ¢f 08 BGEQ  GEd_T0_0.2 : branch if ix! >= 0 
BS ans co. 
3 iki <0. 
71 \¢ : 
3280 8F 50s B1 3071 lz CMPW RO, #*x3280 ® ith 28 
- 3 compare |X! wit tt- 
0 18 9076 315 BGEQ — GEG_TO_2m28 : branch if 1x! >= 200-28 | 
4 17; 
07 18 ; IK! < 2e2-28 
007 19 ; 
007 0 
50 08 70 007 1 MOVD  $*#SD_1.0, RO ; RO/R1 = 1.0 | 
06 0078 ¢ RET 3: return with result = 1.0 
007C 
007C a3 j 
007C 5 3 208-28 =< {Xi < 0.25 
007C 63 
007¢ 7 
007C 8 GEQ_TO_2m28: 
50 50 64 007¢ 9 MULD  _RO,RO ;Get ARG**2 for POLYD. 
AG AF 05 «650 «=675—(O07F 0 POLYD RO, #DCOSHLEN-1, DCOSHTAB | 
0084 1 > RO/R1 = SUM(Ci*Xxeei) 
04 0084 35 RET : return with result in RO | 
085 ; 
085 4; 
085 §$ : 0.25 s< ixi 
Boag 8 3 
085 
085 8 GEQ_T0_0.25: 
FF76 CF «6550 ) 36°71) = 0085 9 CMPD RO, D127 96 2 3; compare :X! with 127*ln2 
1A 14 OBA 40 BGTR  GTR_TRAN_T27_[0G_2 > branch if ix! > 127*Ln2 
08C g 3 | 
8 C 43; 0.25 =< iki =< 127#Ln2 
OBE Se3 
00000000" EF 16 08 46 JSB MTHSDEXP_R6 : RO/R1 = DEXP(X) 
FF79 CF 0 71 4 CMPD =s_- RO,-«D 2 POWER_28.5 > See if the DEXP(-!X!) is significant 
7 14 009 48 BGTR. ONE_TERA ONLY : if not skip computation 
52 98 9 67 099 4 DIVD3 RQ, $*#S0_1.0, R2 : R2/R = DEXP(-!k!) 
0 60 0090 50 ADDD ~— R2,_-RO : RO/R1 = DEXP(X) + DEXP(-x) 
OAO 251 ONE_TERM_ONLY: 
50 0080 BF A2 A 2¢ SuBwW #*x0080, RO : RO/R1 = (DEXP(X) * DEXP(-x))/2 
04 A 3: return with result in RO/R1 
Ab 34 
A6 $3 
A6 256: 127#Ln2 < Ix! 


2 
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A 34 : 
A 
i § GTR_THAN_127_L0G 3: 
FFSD CF 9 71 A 60 CMPD” RO, 0_128_L0G_2 ; Check for possible overflow 
15 146 OOAB 261 BGTR ERROR ; | 
50  FF67 CF 6 A 86 SUBD i) Les 2_H1, RO 3 
52 50 Free CF 68 OOBR 564 MULDS Stee RO, R2 | 
50 gf 62 per 65 SUBD R2, RO ~ 
04 00C 96 RET 
¢C 67 ; 
C 98 3 128*ln2 < ix! 
‘$99 
7E OO'SF 9A O0C 71 ERROR: MOVZBL #MTHSK_FLOOVEMAT, -(SP) ; condition value 
50 01 OF 79 O0C i ASHQ #15, #T, RO ; RO = result = reserved operand -0.0 
OCA 7 : goes to signal mechanism vector 
OCA 74 ; UCHFSL_MCH_RO/R1) so error handler 
OCA 75 3 can modify the result. 
00000000'GF 01 #+=FB OOCA 7 CALLS #1, G*MTHSS$SIGNAL 3 gignet error and use real user's PC 
00D1 7 3; independent of CALL vs JSB 
04 B00) 8 RET ; return = RO/R1 restored from CHFSL_MCH_RO/ 
+44 80 
OD 81 
00D 82 .END 


MTHSDCOSH 


3 
3; Double Floating Hyperbolic Sains rout 16-SEP-19 
Symbol table oo Be 94:3 


SoSEP-1986 1:21:57 EMTARTL.S 


DCOSHLEN = 
DCOSHTAB R 1 
D_127_L0G_ R 1 
D_128 LOG_ R 1 
D2 POWER 28.5 R 1 
D-L6G_2_HT J a 1 
D_LOG_ gf R 1 
ERR R 1 
geQn10~ cR Ol 
GTR “THEN ns, -L06_2 0 AG R it 
L = 00000004 
MTHSSJACKET_HND teeeeere =X 01 

THSSSIGNAL eeeceeee §=X 00 
MTHSOCOSH 00000058 RG 1 
MTHSDEXP_R6 eeeerere 8X 99 
MTHSK_FLOOVEMAT eerreree x 
ONE _TERM_ONLY QO0000A0 R 01 

$0_T.0 = 00004080 
vACUE = 00000004 

eon wee mee erm ene ae be 
: ; Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
ABS 00000000 < 0.) 00 ¢ 0.) NOPIC USR CON ABS LCL 

MTHSCODE 000000D2 ¢ 210.) O1¢ 1.) PIC USR CON REL LCL 


ee ee 


CPU Time 


Page faults 


Initialization + :00:00.10 0:00:00.7 
commend processing vy 700: 3-88 * 0: 1M. $3 
Syabol table sort 0 : : a7 0: ie -09 
Symbol table output : 03 0: 8; -08 
Psect synopsis output 0: 


Cross-reference output 
Assembler run totals 


of BRR RS 


Ibs working + inte was 750 
120 ‘ memory were used to buffer th 


bytes pages) of virtua intermediate code. 
Th re were 10 pages of symbol table space allocated to hold non-local and 
source Ht A were read in Pass 1, producing 11 object records in Pass 2. 
1 page of virtual memory was used to “define 1 macro. 


Local symbols. 


v04-00 Page 
MTHDCOSH.MAR: 1 


— — NORD NOWRT NOVEC BYTE 


RD NOWRT NOVEC LONG 


i 
— 
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VARS Ra Macro Run Statistics MTHRTL sre ATHDCOSH. MAR; 1 (4) | 


teem new m were eer er eco see ce een + 


! Macro Library statistics ! 


ger eer ere eoe nm w eee ewer ecereere = $ 


“$255$DUA28: CSYSLIBJSTARLET.MLB;2 0 
0 GETS were required to define 0 macros. 


| 
Macro library name Macros defined | 


There were no errors, warnings or information messages. 
MACRO/ENABLE=SUPPRESSION/D1 SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHDCOSH/OBJ=OBJ$:MTHDCOSH MSRC$:MTHJACKET/UPDATE=(ENHS :MTHJACKET) #MSRC 
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